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Abstract 

An  imaging  system  that  consists  of  a laser  scanning 
a surface  and  a heterodyne  receiver  that  measures  the  back- 
scattered  field  is  considered.  When  the  scanned  surface  is 
rough  compared  to  the  wavelength  of  the  incident  laser  beam 
the  coherent  properties  of  the  laser  beam  are  destroyed  in 
the  backscattered  field.  This  incoherence  induces  a noise 
in  the  resulting  image  that  is  commonly  referred  to  as 
"speckle". 

The  rough  surface  is  modeled  crudely  by  multiplying 
the  incident  scalar  field  by  a reflectance  term  and  a random 
phase  term.  The  reflectance  is  the  "signal"  that  is  desired 
to  be  measured.  The  two  dimensional  fields  are  propagated 
from  one  plane  to  another  through  the  Huygens-Fresnel  inte- 
gral. The  random  phase  is  considered  to  be  a zero  mean 
stationary  Gaussian  random  process  whose  variance  and  cor- 
relation distance  are  a function  of  the  rough  surface.  It 
is  shown  that  the  correlation  function  of  the  field  is  very 
narrow  when  the  field  is  reflected  from  surfaces  that  are 
rough  compared  to  the  optical  wavelength.  A complex  cur- 
rent representation  is  used  to  show  that  the  mean  of  the 
output  from  the  optical  detector  is  zero.  However,  the 
statistics  of  the  amplitude  and  amplitude  squared  of  the 
current  do  result  in  a mean  "signal"  and  they  are  developed 
in  a manner  similar  to  the  well  known  narrowband  noise 
model.  Second  moment  models  for  detection  of  the  amplitude 
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of  the  current  and  for  detection  of  the  amplitude  squared 
of  the  current  are  presented.  The  mean  "signal"  and  co- 
variance  "noise"  functions  are  related  to  the  field  corre- 
lation function,  the  reflectance  of  the  surface,  and  the 
system  parameters.  The  system  parameters  Include  the 

scanning  velocity,  the  Gaussian  laser  beam  spot  size,  the  | 

receiver  aperture  size,  the  optical  wavelength,  and  the 

observation  distance.  The  noise  models  describe  both  the 

average  speckle  cell  size  In  the  image  and  the  contrast 

in  the  image  that  is  due  to  the  speckle  noise.  It  Is 

shown  that  the  form  of  the  model  Is  the  same  for  both  the 

far  field  and  near  field  cases.  The  models  developed  In 

this  thesis  provide  a basis  for  determining  the  "optimum" 

signal  processing  method  for  producing  the  "best"  Image 

quality. 


A SPECKLE  NOISE  MODEL  FOR  OPTICAL 
HETERODYNE  LINE-SCAfi  IMAGERY 

I , Introduction 

The  objective  of  this  thesis  Is  to  analyze  an  optical 
heterodyne  line-scan  Imaging  system  for  the  case  of  rough 
surface  reflection.  Fig,  1 Is  an  Illustration  of  the  system 
with  an  accompanying  simple  block  diagram.  The  models 
developed  In  this  thesis  will  Include  all  but  the  final  two 
blocks  shown  In  Fig.  1.  In  the  system  a laser  Illuminator 
scans  the  object  surface  In  a pattern  that  will  result  In 
complete  "once  only"  coverage  of  a certain  area.  The 
scanning  system  also  reflects  the  field  backscattered  from 
each  Illuminated  spot  onto  an  optical  heterodyne  detector  as 
the  beam  Is  scanned  over  the  surface.  The  backscattered 
field  Is  a function  of  the  reflectivity  of  the  object's 
surface  therefore,  an  image  of  the  surface  can  be  obtained 
from  the  detected  field.  If  the  object's  surface  Is  opti- 
cally smooth  then  specular  reflection  will  occur  and  the 
field  will  not  be  reflected  directly  back  to  the  detector 
except  when  the  beam  Is  normally  incident  on  the  object  sur- 
face. However  most  surfaces  are  rough  with  respect  to 
optical  wavelengths  so  that  a portion  of  the  field  Is  con- 
stantly backscattered  to  the  detector.  While  this  rough 
surface  provides  the  energy  needed  for  detection  It  also 
distorts  the  reflected  wave  and  the  resulting  Image  because 
the  exact  form  of  the  surface  roughness  Is  "a  priori" 
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Fig.  1.  The  Optical  Heterodyne  Line-scan  Imaging  System 


unknown  and  constantly  changing  as  the  beam  Is  moved  to 
different  areas.  This  image  distortion  is  commonly  called 
"speckle" . 

In  this  thesis  all  optical  fields  will  be  considered 
monochromatic  or  quas imonochromati c and  will  be  mathemati- 
cally propagated  from  one  point  to  another  by  using  the 
Huygens-Fresnel  integral.  The  rough  surface  will  be  modeled 
using  statistics  and  the  optical  detector  output  current 
will  be  based  on  known  detector  models.  The  resulting 
statistics  of  the  detector  output  current  will  be  used  in 
modeling  two  methods  of  signal  detection.  The  statistics 
of  the  final  detected  signal  will  be  a function  of  the 
system  parameters  such  as  the  laser  beam  spot  size,  the 
scanning  velocity,  the  detector  aperture  size,  and  the  opti- 
cal local  oscillator  field.  The  effect  of  the  parameters  on 
the  reflectivity  "signal"  and  the  speckle  "noise"  will  be 
discussed  within  the  context  of  each  model. 

First  the  background  necessary  for  development  of  the 
system  model  will  be  presented.  Then  the  system  model  will 
be  developed  and  finally  the  two  types  of  signal  detection 
models  will  be  presented  and  discussed. 
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II  BACKGROUND  THEORY 


Before  the  system  model  can  be  developed  It  1s 
necessary  to  establish  a background  In  several  areas.  In 
order  of  discussion  thes?  areas  are:  the  phenomenon  of 
"speckle,"  the  forms  of  the  Huygens-Fresnel  Integral,  the 
use  of  statistics  In  modeling,  and  the  process  of  detection. 
Speckle 

When  the  Image  of  an  object  is  produced  either  through 
direct  detection  or  heterodyne  detection,  a random  inter- 
ference pattern  also  results  which  breaks  up  and  distorts 
the  desired  image.  This  random  interference  pattern  has 
been  named  "speckle"  and  is  the  result  of  the  roughness  of 
most  object  surfaces  compared  to  the  wavelength  of  the  inci- 
dent radiation.  Of  course,  the  radiation  must  be  coherent 
to  observe  an  interference  pattern.  The  most  obvious  example 
of  speckle  is  the  direct  observation  of  the  pattern  produced 
at  a rough  surface  that  is  being  illuminated  by  a visible 
laser.  This  phenomenon  was  first  noted  in  the  early  1960's 
by  Oliver  (Ref  1:220)  and  by  Rigden  and  Gordon  (Ref  2:2367- 
2368).  The  use  of  the  term  "speckle"  has  grown  to  include 
any  random  diffraction  pattern  which  disrupts  the  coherent 
properties  of  some  form  of  radiation.  Speckle  occurs  in 
.many  different  circumstances.  Including  rough  surface  scat- 
tering, atmospheric  transmission,  and  holography.  Much 
Information  covering  the  different  areas  of  the  subject  is 
now  available  including  a complete  Journal  of  the  Optical 
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Society  of  America  volume  (Ref  3)  and  a book  edited  by  Dainty 
(Ref  4).  Dainty  also  Includes  a good  Introduction  that 
briefly  covers  the  history  of  speckle  related  phenomena. 

Also  the  theory  of  the  scattering  of  electromagnetic  waves 
from  rough  surfaces  Is  excellently  described  in  a book  by 
Beckman  (Ref  5).  Goodman  has  specialized  much  of  this  Infor- 
mation for  the  case  of  optical  Intensity  radiation  (Ref  4{ 
Chapter  2,  16:1688-1700),  The  major  concern  In  this  thesis, 
however,  will  be  for  the  case  of  the  direct  measurement  of 
optical  fields. 

Huygens-Fresnel  Integral 

The  real  scalar  optical  field,  u(f,t),  Is  In  general  a 
function  of  time,  t,  and  space,  r,  where  r Is  a vector  with 
Cartesion  components  (x,y,z).  A temporally  modulated  field 
located  about  frequency,  f^,  can  be  written  as 

u(r,t)  = A(r,t)cosl2iif^t  + <})(?, t)] 

" Re{A(r,t)expI-j4>(r,t)]expI-J2nf^t]} 

» Re{U(r,t)exp[-j2nfQt] } (1) 

where  Re{*}  is  the  real  operator  and  U(r,t)  is  the  complex 
envelope  of  the  complex  representation  of  u(F,t).  The  tem- 
poral Fourier  transform  of  u(r,t)  Is 

ft[u(r,t)J  “ ^F^IU(r,t)e’J^’'%b 

+ jFtlU*(r,t)eJ2’'^o^]  (2) 

00 

where  F^Ig(t)]  « /g(t)expl -J2TTf  tldt  and  the  asterisk,  (*), 


represents  the  conjugation  process.  From  a graphical  repre 
sentatlon  of  Eq.  (2),  shown  In  Fig.  2,  It  can  be  seen  that 
the  transforms  of  U(r,t)  and  U*(f,t)  are  centered  at  the 
minus  and  plus  frequencies  of  f^  respectively.  If  u(r,t) 


Is  unmodulated  the  complex  envelope  equals  U(r)  and  the 
Fourier  transform  becomes  two  Dirac  delta  functions  that 
are  located  only  at  frequencies  plus  and  minus  f^.  This 
field  Is  then  called  monochromatic  because  of  the  one 
frequency  component.  If  u(r,t)  Is  modulated,  I.e.,  time 
varying  amplitude  or  phase  as  represented  In  Fig.  2,  then 
the  field  Is  called  quasimonochromatl c provided  the  band- 
width, B,  of  F.[U(r,t)]  Is  much  less  than  f , I.e.,  B<<f  . 

w U O 

Thus  for  either  monochromatic  or  quasimonochromatl c fields 
the  complex  envelope  of  u(r,t)  can  be  used  to  represent  the 
field  and  the  real  field  can  be  obtained  through  Eq.  (1). 
Therefore,  throughout  most  of  this  thesis,  the  complex  en» 
velope  will  be  used  and  the  expont^'ntlal  optical  time  func- 
tion suppressed.  In  addition,  the  field  Is  often  described 
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at  a particular  z plane  so  that  the  field  varies  only  In 
X and  y.  Thus,  the  unmodulated  optical  time  varying  field 
at  a z plane  can  be  described  through  the  complex  envelope, 
U(x.y). 

As  shown  by  Goodman  (Ref  7:58),  the  complex  envelope  of 
a monochromatic  field  at  a point  in  one  plane  can  be 
expressed  in  terms  of  the  complex  envelope  of  a field  at 
another  plane  through  the  following  Huygens-Fresnel  integral 
equation: 


U(x,y)  » // 


z 


expUkr^i  ]cos(n,rn,  )U(a,6)dade 


The  vector  is  shown  in  Fig-  3,  k is  the  magnitude  of 
the  propagation  vector  which  is  equal  to  2n/X,  X is  the 
optical  wavelength,  and  cos  is  the  cosine  of  the 

angle  between  Fq^  and  the  normal  to  the  a, 6 plane.  Eq.  (3) 
can  be  simplified  through  several  approximations.  If  the 
angle  between  n and  Fq^  is  limited  to  less  than  18°  then 
the  cos  (n,FQj)  is  approximately  equal  to  one.  The  inte- 
gral over  the  aperture,  Z,  can  be  replaced  by  an  integral 
of  infinite  limits  if  the  finite  extent  of  the  field, 
U(a,8},  due  the  aperture  is  included  in  the  mathematical 
description  of  U(a,p).  In  addition,  if  the  distance  z is 
much  greater  than  the  maximum  linear  distance  from  the  z 
axis,  then  in  the  denominator  can  be  approximated  by  z. 
The  above  restrictions  have  put  the  Integral  in  a simpler 
form  and  are  discussed  in  greater  detail  by  Goodman, 
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Fig.  3.  Diffraction  Geometry  (Ref  7: ^ 


(Ref  7:58),  The  integral  now  becomes:  j 

CO  ^ 

U(x,y)  « //  U(a,e)exp[JkrQ^]dadp  (4)  j 

The  integral  is  still  not  very  useful  because  of  the  rQj 
term  in  the  exponential.  But,  because  the  exponential  is 
very  sensitive  to  small  changes  in  the  value  of  exponent, 
extra  care  must  be  used  in  making  any  approximations  in  the 
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exponential.  The  two  approximations  commonly  made  at  this 
point  are  th.e  Fresnel  and  Fraunhofer  approximations. 


Fresnel  Approximation.  The  Fresnel  approximation  is 
made  by  writing  rp^  from  Fig.  1 as: 

rgj  - [z^  + (a-x)^  + (3-y)^]^^^*z[l  + (5^)^  + (^)^]^/2  (5) 

Since  the  absolute  value  of  the  difference  terms  divided  by 
z is  less  than  one,  the  binomial  expansion  to  the  square 
root  can  now  be  used  to  write: 

*'01  = ^ + 1 (6) 

Now  if  the  maximum  linear  distance  is  such  that 

z'  » ^ [(a-x)2  . (7) 


then  the  higher  order  terms  beyond  the  squared  terms  may 
be  neglected  and  Eq.  (4)  reduces  to  the  Fresnel  approxima- 
tion given  as: 


U(x,y)  = u(a,0)exp[^[(a-x)2 


+ {e-y)'^]]dad8 


(8) 


or 

U(x,y) 


exp[ jkz]exp[|j(x^+y^) 
jXz 

exp[^^^||^  (xa+y6)ldad6 


ff  U(a,6)[exp  -^(a^  + B^)] 

»00 

(9) 


Fraunhofer  Approximation.  The  Fraunhofer  approximation 
Is  used  when: 
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(10) 


max 


If  Eq.  (10)  holds,  then  the  squared  phase  terms  may  be  neg- 
lected In  the  Integral  of  Eq.  (9)  and  the  Fraunhofer 
approximation  is  written  as: 


U(x.y)  = 


exp[jkz]exp[J-y(x^+y^)  ] ” 

= r r 


jXz 

exp[^||^  (xa+  yB)]dadB 


I!  U(a.6) 


(11) 


The  integral  in  Eq.  (11)  can  now  be  recognized  as  the  two 
dimensional  Fourier  transform  of  U(a,B)  evaluted  at  the 
spatial  frequencies  ~ Xz  ^y  " X^‘ 

One-Dimensional  Huygens-Fresnel  Integral.  In  many 
cases  the  effect  of  the  basic  parameters  in  a system  can  be 
adequately  described  in  simpler  terms  through  the  use  of  the 
one-dimensional  Huygens-Fresnel  intergral.  In  this  case  the 
field  at  the  z plane  is  variable  in  only  one  lateral  direc- 
tion. The  integral  is  now  written  as  (Ref  22:316): 


U(x) 


exp[-J2nJ^ld«  • (12) 

Up  until  this  point  it  has  been  implied  that  the 
Huygens-Fresnel  integral  can  be  used  only  for  monochromatic 
waves.  However,  the  integral  may  be  used  for  quas imonochro- 

I t V*  I 

matic  waves  provided  g >>  ' Ol'max  where  c is  the  speed  of 
light  (Ref  7:108).  B is  the  bandwidth  of  the  complex  en- 
velope as  shown  in  Fig.  2. 
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Propagation  and  Reflection  of  a Laser  Beam.  When 
Eq.  (11)  Is  used  to  propagate  the. output  from  a laser 
source  the  result  Is  a Gaussian  spherical  wave  (Ref  8:306) 
of  the  following  form 

2 2 

U(o,3)  = A exp(  jkz]exp(  (a^  + 6^  ) ]expl^^-4^ — ] (13) 

where  A Is  an  amplitude  term  . The  beam  spot  size  at  the  z 
plane,  w(z).  Is  the  radial  distance  where  the  total  beam 
amplitude  Is  e"^  times  the  center  amplitude.  If  It  is  now 
desired  to  reflect  this  wave  directly  back  from  normal 
Incidence  on  a smooth  surface  using  the  same  coordinate 
system,  then  the  exponential  phase  terms  must  be  conjugated. 
By  conjugating  the  expfjkz]  term,  the  wave  now  propagates 
in  the  oposite  direction,  and  by  conjugating  the  explj^ 

(a  +3  )]  term,  the  spherical  phase  fronts  are  inverted. 

The  Inversion  of  the  wavefront  can  be  seen  by  simply  drawing 
rays  and  th.en  reconstructing  the  wavefront.  More  formally 
the  spherical  wave  can  be  broken  up  into  an  infinite  sum  of 
plane  waves  of  different  angles  of  incidence.  The  surface 
boundary  conditions  can  be  applied  to  each  plane  wave  and 
the  resulting  sum  then  used  to  reform  the  reflected  spheri- 
cal wave.  In  addition,  since  a wave  Is  seldom  totally 
reflected,  a reflection  coefficient,  a (a, 8),  that  Is  a 
•real  function  of  the  particular  surface  must  be  multiplied 
times  the  incident  field  to  complete  the  process  of  reflec- 
tion, (Ref  9:74).  In  general  the  reflection  coefficient 
depends  on  the  Incident  angle  and  polarization  of  the 
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Incident  field  and  on  the  index  of  refraction  of  the  surface. 
The  reflected  field  is  modeled  as 

U^(a.0)  = a (a,B)U*(a,e)  (14) 

where  r denotes  the  reflected  wave  and  i denotes  the  inci- 
dent wave 

Statistics  in  Modeling 

In  many  instances  the  exact  nature  of  the  desired  sig- 
nal field  is  unknown  either  because  the  signal  itself  is 
unknown  or  because  a known  signal  has  been  distorted  by 
some  type  of  unknown  interference,  or  both.  In  these  cases 
it  is  useful  to  model  the  system  through  the  use  of 
statistics. 

In  statistics  a random  variable  is  an  unknown  function, 
X,  whose  possible  values  are  best  described  by  a probabil- 
ity density  function  (Ref  10:92-136).  Each  possible  value 
of  X is  called  a sample  point.  In  many  cases  x is  also  a 
function  of  time,  termperature , space,  or  any  other  index 
and  is  then  called  a random  process  (Ref  10:298-339).  Thus 
a random  process  is  a collection  of  random  variables  that 
are  indexed  by  another  parameter  so  the  probability  density 
also  becomes  a function  of  the  index.  For  any  particular 
value  of  the  index,  (x)  is  a random  variable.  For  any  fixed 
value  (sample)  of  a random  variable  (x)  is  still  a function 
of  the  index  and  is  called  a sample  function.  To  describe 
the  random  process  completely,  all  orders  of  the  combined 
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probability  densities  for  all  possible  values  of  the  index 
must  be  known  (Ref  10:311).  But,  in  practice  this  requires 
more  knowledge  of  the  random  process  than  is  usually  avail- 
able. Therefore  a limited  statistical  description  of  the 
system,  called  a second  moment  model,  is  commonly  used. 

The  mean,  variance,  correlation,  and  covariance  function 
are  used  in  a second  moment  analysis. 

The  mean  of  a random  process,  x,  is  in  general  a func- 
tion of  the  index,  t,  and  is  denoted  as  where  El*] 

is  the  expected  value  operator  (Ref  lOrChapter  7).  The 
mean  is  the  first  moment  of  the  process  and  is  often  called 
the  expected  value  of  the  process  or  the  ensemble  average. 
The  mean  crudely  represents  the  "most  likely"  sample  func- 
tion of  the  random  process  and  is  usually  directly  related 
to  the  desired  signal.  The  variance  describes  the  average 
of  the  square  of  the  fluctuations  about  the  mean  and  is 
defined  by  (Ref  10:244) 

a^(t)  = E[(x^-E[x^.])^]  = E[x^]  - E^(x^]  (15) 

where  again  the  index,  t,  has  been  included  because  x is  a 
random  process.  The  variance  is  usually  considered  an 
Indication  of  "noise"  in  the  system  because  it  is  a measure 
of  the  variation  or  fluctuations  of  the  random  process  from 
the  desired  mean  term.  The  square  root  of  the  variance  is 
usually  denoted  as  the  standard  deviation  or  the  rms  varia- 
tion from  the  mean. 
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since  the  process  varies  with  the  Index,  the  mean  and 
variance  are  not  sufficient  Information  to  describe  how  the 
process  is  related  to  Itself  at  different  values  of  the 
index.  This  relationship  Is  described  to  some  degree 
through  the  correlation  Cor  autocorrelation)  and  covariance 
functions.  The  covariance  function  Is  defined  as  (Ref  10: 
317) 


CxCti.tg)  = E[(x^^-E[x^^])(x^^-EIx^^]) 


® EIx^  X.  ] - E[x.  ]E[x^.  ] 

ti  tg  tj  t2 


where  called  the  correlation  function.  If 

^1~^2  covariance  is  equal  to  the  variance  and  if 

the  mean  is  zero,  the  covariance  is  equal  to  the  correla- 
tion function. 


Because  the  covariance  depends  on  the  individual 
values  of  the  Index,  it  may  still  require  more  information 
than  Is  usually  known.  But,  In  some  cases  the  statistics 
can  be  considered  stationary  at  least  at  the  second  moment 
level  (Ref  10:325).  Under  this  condition  the  mean  becomes 
a constant  (Independent  of  the  Index)  and  the  covariance 
and  correlation  functions  depend  only  on  the  absolute 
•difference  In  the  Index  values,  t2-t^>»At«T,  This  assumption 
may  appear  Invalid  for  a particular  random  process,  but  It 
only  means  that  the  process  Is  being  described  by  a rather 
crude  but  still  useful  set  of  statistics.  Under  the  condl- 
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tion  of  statlonarity,  Eq.  (16)  now  becomes; 

c^Ct)  « R^(t)  - E^Ix]  ■ a7) 

The  covariance  and  correlation  functions  of  a real  random 
process  can  be  shown  to  be  even  functions  of  t with  their 
maximum  magnitude  at  t=0  (Ref  10:323),  Also,  for  a nonperi- 
odic process,  the  correlation  function  approaches  the  square 
of  the  mean  and  the  covariance  approaches  zero  as  x 
approaches  Infinity,  In  many  cases  the  effective  non-zero 
value  of  the  covariance  occurs  over  just  a short  length  of 
T from  the  origin.  This  period  of  x Is  called  the  correla- 
tion or  coherence  period  (time,  distance),  x^,  of  the  random 
process.  Although  the  exact  point  at  which  the  covariance 
can  be  considered  to  be  zero  for  a particular  function  may 
be  rather  arbitrary,  x^  roughly  represents  the  distance  (or 
time  separation)  between  samples  of  the  process  such  that 
the  samples  become  unrelated  I.e,,  uncorrelated,  A typical 
sample  function  with  statistical  parameters  crudely  Ident- 
ified Is  shown  In  Fig.  4a.  A sample  correlation  function 
is  shown  In  Fig.  4b. 

Measurement  of  Optical  Signals 

To  retrieve  the  Information  from  an  optical  signal, 
there  must  be  some  means  of  measurement.  Most  available 
optical  detectors  are  square  law  devices  that  measure  the  i 

intensity  of  the  Incident  field,  I.e.  |u(x,y)f^  (Ref  11:1819), 

However,  homodyne  and  heterodyne  techniques  can  be  used  to 
measure  the  actual  field  with  Its  accompanying  phase 
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Information.  While  the  system  to  be  discussed  in  this 
thesis  is  a heterodyne  system,  a brief  discussion  of  direct 
detection  will  be  given  first  so  that  the  difference  In 
the  necessary  field  statistics  can  be  related.  The  detectors 
considered  here  will  be  modeled  as  ideal,  i.e.,  all  quantum 
effects  and  noise  terms  will  be  neglected.  This  will  allow 
the  speckle  noise  term  to  be  isolated  in  this  thesis.  The 
detector  noise  terms  can  always  be  added  to  the  model  if  it 
becomes  necessary. 

Direct  Detection.  The  ideal  instantaneous  output  of 
an  optical  detector  is  defined  as  CRef  12:91) 

Kt)  = ^ //  |U5(x,y)|2dxdy  (18) 

0 Aj 

where  q is  the  charge  of  an  electron 

n Is  the  quantum  efficiency  of  the  detector 

J 

h is  the  Planck's  constant 

f^  Is  the  optical  frequency 

P^  is  the  optical  power  collected  by  the  detector 
surface 

Ajj  Is  the  area  of  the  detector  surface 

UjCx,y)  is  the  signal  field  incident  on  the  detector. 

Thus  this  sytem  directly  measures  the  intensity  of  the 
field  and  Is  referred  to  as  direct  detection. 

The  mean  current  is 

ElUt)]  = ^ //  EIUjx,y)U*(x,y)]dxdy  (19) 

^d 
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where  the  asterisk  denotes  the  complex  conjugate  of  the 
field.  Therefore  the  expected  value  of  the  current  Is  re- 
lated to  the  expected  value  of  the  intensity  of  the  field 
or  the  correlation  function  of  the  field.  Similarly  the 
correlation  function  of  the  current  is  related  to  a fourth 
order  moment  of  the  field  as  shown  below; 

EIi(t)i{t')]  = ////  ElU5(x,y)U*(x,y)U^(x',y')U* 

^d 

(x' ,y’ )]dxdydx' dy‘  (20) 

Heterodyne  Detection.  In  a heterodyne  system  the  in- 
coming signal  field  is  added  to  a local  oscillator  field, 
means  of  a beam  splitter.  The  local  oscil- 
lator field  is  at  an  optical  frequency  of  With 

thi'  arrangement  the  ideal  detector  output  current  becomes 
(Ref  13:481-487) 

1(t)  = lUg(x,y,t)  + ULQ(x,y,t)  l^dxdy 

^d 

= ITfJ  { lUgCx.y)  l^+2Re[U^Cx,y)u[Q(x,y)exp[-j2nf  jpt]] 
^d 

+ lU^QCx.y) l^ldxdy  (21) 

where  Re[*J  Is  the  real  operator.  The  first  and  last  terms 
In  the  Integral  are  centered  at  zero  frequency  and  can  be 
.filtered  out  so  that  the  remaining  output  is  at  the  inter- 
mediate frequency,  f^p,  as  shown  below 

^ Re{Ug(x,y)U*Q(x.y)exp[-J2nfjpt]}dxdy  (22) 
^d 
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since  the  expected  value  and  real  operator  are  both  linear 
they  can  be  Interchanged  and,  assuming  that  the  signal  field 
Is  the  only  random  quantity,  the  expected  value  of  the  cur- 
rent Is: 

E[1(t)]  = //  Re{E[Ug(x,y)JU*p(x,y)exp[-2jnfjptJ}dxdy  C23) 

^d 

Thus  for  heterodyne  detection  the  mean  of  the  current  is 
related  to  the  mean  of  the  field.  This  is  in  contrast  to 
the  case  of  direct  detection  where  the  mean  of  the  current 
from  Eq.  (19)  is  proportional  to  the  correlation  of  the 
field,  i.e.,  E[Ug(x,y)U*(x,y)] . The  correlation  function 
of  the  current  for  heterodyne  detection  is: 

2 

EM(t)i(t')]  = (^)  ElfffJ-  Re{Uj(x,y)U*Q(x,y) 

^d 

expl-j2nfjpt  }RelU^(x' ,y' IU^qCx' ,y' ) 

expl-  j2nf  j pt ' ] dxdydx ' dy ' 

2 

" Re{////  EIUj(x,y)U5(x',y')]U*Q(x,y) 

I 

Ulq(x' ,y' )expI-j2nfjp(t+t')]+ElU3{x,y)  | 

U*Cx'  ,y'  )JULQ(x,y)U^^Q(x'  ,y'  )exp[-J2nf  jp 
(t-t ' ) Jdxdydx ' dy ' } C24) 

where  identity  RelAlRelB]  RelAB+AB*]  has  been  used. 

Therefore  in  heterodyne  detection  the  correlation  of  the 
current  Is  related  to  the  correlation  of  the  real  parts  of 
the  field  and  is  in  general  nonstationary. 
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The  discussion  of  the  background  theory  Is  now  com- 
plete. The  task  is  to  now  use  these  models  and  ideas  that 
have  been  presented  in  this  section  to  develop  a complete 
system  model  for  the  heterodyne  line-scan  Imaging  system. 
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III.  THE  SYSTEM  MODEL 
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Now  that  the  background  theory  has  been  discussed  the 
system  model  can  be  developed.  First  a phase  model  for  the 
rough  surface  will  be  established,  then  this  phase  model 
will  be  incorporated  into  the  total  heterodyne  line-scan 
imagery  model.  From  the  system  model  the  mean  and  covar- 
iance functions  of  the  complex  output  current  will  be 
determined.  For  simplicity  the  model  will  be  developed  in 
one  dimension  as  explained  later. 

The  Rough  Surface  Phase  Model 

The  reflection  ^rom  a rough  surface  is  in  general  a 
very  complicated  process  that  involves  the  surface  reflec- 
tion coefficients,  the  surface  height  variations,  the 
macroscopic  and  microscopic  angles  of  incidence,  and  the 
polarization  of  the  incident  field.  Since  the  system  dis- 
cussed here  involves  the  imaging  of  many  different  types  of 
surfaces  with  characteristics  beyond  the  system's  control, 
and  since  the  usable  information  is  only  in  the  direct  back- 
scattered  radiation,  a more  simple  phenomelogical  model  will 
be  used.  The  model,  which  was  also  used  in  a similar  direct 
detection  imaging  problem  (Ref  14:779-785),  is  less  com- 
plicated than  Beckman's  (Ref  5:Chapter  5)  but  is  closely 
related.  The  model  will  be  developed  under  the  following 
conditions:  (1)  scalar  fields  will  be  used,  (2)  depolari- 

zation effects,  multiple  scattering,  and  shadowing  will  be 
neglected,  as  is  typically  the  case  (Ref  5:Chapter  3 and  5, 
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6:1689)  since  the  effects  are  often  small  and  very  difficult 
to  describe  mathematically,  (3)  the  surface  will  be  con- 
sidered rough  as  compared  to  the  optical  wavelength  as  dis- 
cussed In  more  detail  later  so  that  backscattered  radiation 
exists  and  total  specular  reflection  does  not  occur,  (4) 
the  radius  of  curvature  of  the  surface  Irregularities  will 
be  assumed  large  as  compared  to  the  optical  wavelength 
(Ref  5:20),  and  (5)  the  field  reflectivity  will  be  considered 
a real  function  of  space  with  a value  between  zero  and  one 
so  that  it  attenuates  the  Incident  field  such  that  the 
reflectivity  is  the  surface  characteristic  to  be  measured, 

1 .e.  the  signal . 

Under  the  above  conditions  the  reflection  from  a rough 
surface  will  be  modeled  by  including  a random  phase  term  in 
the  process  of  reflection  from  a smooth  surface  given  by 
Eq.  (14)  In  the  background  section.  The  reflected  wave  1s 
then : 

U^(a)  = a(a)e^®^'‘^uj(a)  (25) 

The  random  phase  can  now  be  related  to  the  surface  heights. 
First,  the  surface  heights,  h(a),  will  be  modeled  as  a sta- 
tionary zero  mean  Gaussian  random  process.  The  zero  mean 
follows  from  the  fact  that  any  constant  reference  can  be 
added  to  the  model  such  that  the  ensemble  average  becomes 
zero,  A Gaussian  distribution  Is  commonly  used  to  describe 
a rough  surface  (Ref  4:65,  3:1153,  1195,  1205,  1212,  1224, 
5:80)  but  may  break  down  for  some  polished  man  made  surfaces. 
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As  win  be  discussed  later  tMs  Is  not  a particularly 
! necessary  assumption,  but  ft  does  make  the  problem  easier 

! mathematically, 

I 

1 

I As  can  be  seen  from  the  sample  function  of  Ftg.  5,  the 

total  path  difference  between  a wave  reflected  at  hCa)  and 
a wave  reflected  at  hCo)  referenced  to  a constant  spherical 

i * 

j phase  front  Is 

1 2 2 ■ 2 2 

Ar  “ d+^22  ” h(a)+d— " h Ctx ) "d d 2 2 ° - 2 h Co ) (26) 

an  so  the  phase  variation  becomes; 

0(a)  « ^(2h(a))  = ^ h(al  (27) 


Induced  by  the  Rough  Surface 


( 

J 
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The  mean  of  the  phase  Is 


EieCal]  = ^ EIhCol]»0 


C28i 


and  the  variance  of  the  phase  is: 

9 ATT  ^ 9 ^ 

cl  ■ ci|)  cl  ■ 


C29} 


From  Eq.  (29)  it  can  be  seen  that  the  rms  roughness  of  the 
surface,  can  be  directly  compared  to  the  wavelength, 
When  the  wavelength  is  large  compared  to  a^,  the  rms  varia- 
tion  of  the  phase  is  small  and  there  would  be  little  pro- 
blem with  interference.  But,  as  approaches  the  value  of 
X or  greater,  the  rms  phase  variations  become  larger  than 
2II  radians,  and  destructive  interference  will  result.  An 
example  of  some  values  of  rms  surface  roughness  is  shown  in 
Table  1.  Most  exterior  metal  surfaces  will  fall  in  the 


Table  I Values  of  Surface  Roughness 


Name 

rms  roughness 

average  peak  to 
peak  height 

micro- 

inches 

microns 

micro- 

inches 

microns 

Rough 

1000 

25.4 

3500 

88.9 

Semi  rough 

500 

12.7 

1750 

44.45 

Medium 

250 

6.35 

875 

22.2 

Semi  fine 

125 

3.18 

455 

11.56 

(From  Ref.  15:859) 
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category  of  rough  to  semlrough  and  thus  would  also  be  con- 
sidered rough  as  compared  to  an  optical  wavelength.  Also 
most  other  nonmetallc  surfaces  are  rough  compared  to  metal 
surfaces  so  they,  too,  would  be  rough  compared  to  an  opti- 


cal wavelength. 


Later  In  determining  the  moments  of  the  fields  It  will 


be  necessary  to  take  the  expected  value  of  the  random 
exponential  phase  terms  as  follows: 


£ I e J ^ ® ^ ^ ^ ] a Eje^"^! 


The  above  two  equations  represent  a form  of  the  character- 
istic function  (Ref  10:419) 


(|)^(v)  =»  Ele^^’'] 


for  the  case  of  v=l.  For  a Gaussian  random  variable,  x. 


the  characteristic  function  1s  well  known  and  equals  (Ref 
10:420): 


♦^(v)  - expIjvElxJ 


2 2 


Thus,  for  Eqs  (30),  (31)  and  C32),  v=l  and  the  mean  Is  zero 
so  the  characteristic  function  equals: 

Ele^'*']  e^  (35) 


where  t=0,<|),  or  v as  required.  The  variance  of  G was 
determined  In  Eq.  (29).  For  a stationary  zero  mean 


Gaussian  random  process  the  variance  of  (|>  is  computed  as 
fol 1 ows : 


1 


0^  = E[(6(a)-6(a' ))^]-E^I0(a)-0{a' )] 

= E[d^(a)  +0^(a' )-20(a)0(a' ) ] 

= ZOq  - 2Rg(Aa) 

= 2og(l-p(Aa))  (36) 

where  Aa=o-a'  and  p(Aa)  is  the  correlation  function  of  the 
phase  normalized  to  one  at  Aa=0.  p(Aa)  is  also  identically 
the  normalized  correlation  function  of  the  surface  heights, 
h(a).  Similarly 

= 2a^(l+p(Aa))  (37) 

The  Heterodyne  Detector  Current  Model 

The  rough  surface  model  is  now  complete  and  must  be 
incorporated  into  the  system  model.  The  output  field  from 
the  illuminating  laser  will  be  propagated  to  a rough  sur- 
face, reflected,  and  propagated  back  to  the  detector  by 
means  of  the  Huygens-Fresnel  integral  equation.  The 
detector  current  will  be  determined  from  this. field.  Then 
the  statistical  moments  of  the  current  will  be  determined 
and  it  will  be  shown  that  the  system  can  be  described  in 
terms  of  the  complex  current  output  of  the  detector. 

The  Real  Current.  The  mechanics  of  the  heterodyne  line- 
scan  system  scan  the  laser  radiation  across  the  scene  to  be 
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Imaged  and  the  detector  measures  the  backscattered  field. 

The  system  mechanics  also  keeps  the  detector  surface  and 
local  oscillator  field  aligned  normal  to  the  direction  the 
laser  is  pointing.  However,  the  angle  of  incidence  between 
the  laser  field  and  the  object  surface  will  change  as  the 
beam  is  swept  across  the  object.  In  principle  the  coordin- 
ate rotations  required  to  describe  this  scanning  system 
could  all  be  included  in  the  Huygens-Fresnel  integral  but 
the  problem  becomes  more  complicated  than  useful.  There- 
fore, the  system  will  be  modeled  with  the  laser,  detector, 
and  associated  fields  in  a fixed  coordinate  system.  The 
laser  beam  will  be  normally  incident  upon  an  object  surface 
whose  reflectivity  and  random  phase  characteristics  are 
effectively  moving  beneath  the  beam  in  time  and  with  veloc- 
ity, V.  Both  the  spreading  of  the  beam  spot  and  the  changing 
surface  velocity  that  occur  as  the  beam  is  swept  over  the 
surface  wilj  be  ignored.  But,  the  model  will  be  a function 
of  the  beam  width  and  the  velocity  so  the  above  effects  can 
be  determined  simply  by  varying  these  parameters.  Also,  at 
any  instant  of  time,  the  detector  output  is  proportional  to 
an  average  of  all  the  reflectance  values  at  each  point  with- 
in the  beam  spot  and  detector  field  of  view.  Because  the 
laser  beam  is  circular  in  shape,  the  reflectance  points 
• covered  by  the  width  of  the  beam  at  each  point  in  the 
scanning  direction  will  change  slightly  as  the  beam  is  moved. 
But,  the  basic  change  in  the  average  reflectance  is  caused 
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by  the  reflectance  points  that  enter  and  leave  the  beam  In 
the  scanning  direction.  Thus  the  model  will  be  developed  In 
the  scanning  direction  only,  where  the  reflectance  at  each 
point  In  this  direction  can  be  thought  of  as  an  average  of 
the  reflectance  over  the  width  of  the  beam.  The  extension 
of  the  model  to  both  lateral  dimensions  Is  straightforward 
but  more  tedious. 

The  system  geometry  Is  shown  In  Fig.  6.  From  Eqs  (12) 


Fig.  6 System  Geometry 


and  (13)  the  one  dimensional  field  from  the  laser  at  the  ob- 
ject surface  Is: 
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U^(a)  = 


exp[ j (kz-J) ] 


I^A  expt-(^)'].xplj!^) 


(XZ) 


(38) 


From  Eq.  (25)  where  the  reflectivity  and  random  phase  are 
now  properties  of  a moving  surface,  the  backscattered  field 
is: 


.0(a-vt)  * 

U^(a)  = a(a-vt)e'^  U^(a) 


(39) 


The  reflected  field  is  now  propagated  back  to  the  detector 
through  the  Huygens-Fresnel  integral  of  Eq.  (12)  but  since 
propagation  is  now  in  the  negative  z direction  the  phase 
terms  of  the  Huygens-Fresnel  integral  of  Eq.  (12)  must  also 


be  conjugated  so  that  the  field  at  the  detector  is 

2 

2 


_ expt-j(kz-f)]expl-Jkf^l  . 


U^(a)  = 


(kz) 


1/2 


f U^(a)exp[-j^] 


exp[j2n||]da 


(40) 


From  substituting  Eqs.  (38,  39,  and  40)  into  Eq.  (22)  the 
detector  output  current  becomes 


00 

1p(t)  = Re{exp[-j2(kz-J)]  //  Pp( x ) a^ (a+vt ) 

r g 

‘w(z) 


2 2 

expI(:77TT)^]expIje(a+vt)]exp[-j^]exp[-j^] 


exp[j2n^]U*Q(x)exp  t-j2nf  jptjdxda) 


(41) 


where  Prj(x)  is  the  limiting  detector  aperture  function.  A 


simplication  to  Eq.  (41)  can  be  made  by  realizing  that 
Pq(x)  can  represent  a lens  aperture  and  transfer  function 


T 
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or  just  the  detector  aperture  size.  If  the  local  oscil- 
lator field  Is  a plane  wave  the  x Integral  iii  the  above 
equation  yields  the  Fresnel  diffraction  pattern  of  P[j(x). 
However,  this  can  be  simplified  by  realizing  that  In  most 
cases  Pq(x)  is  narrow  enough  to  satisfy  the  Fraunhofer  type 

approximation  given  by  Eq.  (10).  If  this  is  not  the  case  the 
kx2 

exp[-j— phase  term  can  still  be  negated  either  by  a con- 
jugate lens  transfer  function  (converging  lens)  or  by  an 

Identical  phase  term  from  the  local  oscillator  field.  In 

kx2 

any  case  U|^q(x),  exp [ p j-^^]' and  any  phase  term  in  Pq(x) 
can  be  all  combined  to  equal  one,  so  that  Pq(x)  is  now  just 
the  limiting  aperture  size.  The  x integral  is  now  the 
spatial  Fourier  transform  of  P[j(x)  evaluated  at  f^  = 

The  transform  of  Pq(x)  is  denoted  as: 

Equation  (41)  now  reduces  to 

CO 

i^(t)  * BRe{exp[-j2( kz-J) ]exp[-j2nf jpt]  / h(a) 

-*00 

expl  je(a+vt)  ] a(a+vt)expt-j^^|-]da}  (43) 

where  B = 
wnere  b 

2 

and  h(a)  = PoF^“^®’‘P ^ 

h{a)  will  be  referred  to  as  the  system  function  and  repre- 
sents the  combined  result  of  the  beam  spot  size  and  the 
detector  field  of  view. 
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The  Complex  Current.  As  was  seen  in  Eqs.  C23)  and  (24) 
the  moments  of  the  current  could  be  obtained  by  taking  the 
real  part  of  the  moments  of  the  complex  current.  Thus,  the 
moments  of  the  complex  current  will  now  be  determined.  In 
determining  t..«!se  moments  the  phase  terms  in  front  of  the 
integral  in  Eq.  (43)  can  be  dropped.  It  will  be  shown 
later  that  these  phase  terms  are  indeed  unrelated  to  the 
problem  solution.  Thus  the  useful  complex  current  is  de- 
fined as: 

" 1,  2 
i(t)  = B / h(oi)a(a+vt)expIje(a+vt)lexpl-J-=~]da  (45) 

•.00  ^ 

Since  the  rough  surface  phase  term  is  the  only  random 
term  in  Eq.  (45)  the  mean  of  the  current  using  Eq.  (35)  is 

a 2 « 2 

Eli(t)]  = B exp["^  ;h(a)a(c(+vt)expl-j^]da  (46) 

••00 

/®e)2 

where  exp [-  is  the  value  of  the  characteristic  func- 

tion from  Eq.  (35).  The  value  of  the  characteristic  func- 
tion  for  a surface  with  ^ A in  Eq.  (29)  is  5x10'”  or 

p 

less.  For  typical  values  of  A=.l,  n=l»  and  fA=c=3xlO 
Eq.  (46)  becomes  (5x10"^^)  (1.6x10^)  (L/z)  where  L is  the 
value  of  the  integral  in  Eq.  (46).  Since  h(a),  a(a+vt)  and 
expi'Jy”  ] all  are  terms  with  values  less  than  one,  L Is 
equal  to  the  beam  width  or  less  and  thus  the  mean  current 
is  effectively  equal  to  zero. 
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with  a zero  mean,  the  covariance  of  the  current  equals 
the  correlation  of  the  current.  From  Eq  (24)  the  correla- 
tion involves  the  moments:  EIi(t)i(t')]  and  E [ i ( t ) i * ( t ' ) ] . 
Using  Eq.  (45)  the  first  moment  is: 

00 

E[i(t)i(t')]  = If  h(a)h(a ' )a (a-vt )a (a ' -vt ' ) E [exp{ j [6 

• 00 

{a-vt)+0(a'-vt' )] }]exp[-jY(a^+a'^) ] 

dada'  (47) 

From  Eqs.  (33,  34,  35,  and  37)  it  is  seen  that  the  above 
expected  value  term  is  another  characteristic  function  and 
is  equal  to: 

4>^(1)  = exp[-Og(l  + p(Aa+vAt))]  (48) 

where  (a-vt)-(a' -vt' )-Act-vAt.  The  normalized  covariance 
function  of  the  phase,  p(Aa+vAt),  varies  between  one  and 
zero  so  that  Eqs.  (47)  and  (48)  equal  zero  under  the  same  con- 
ditions as  for  the  zero  mean  current  of  Eq.  (46).  The 
other  moment  is: 

00 

E(i(t)i*(t*  )]  = ll  h(a)h(a' )a(ct+vt)a(a’+vt’ )E[exp{j  [0 

(a+vt)-0(a'+vt' )]}]exp[-jj(a^-a'^)] 

dada'  (49) 

From  Eqs.  (35)  and  (36)  the  above  characteristic  function  is 
equal  to: 
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1' 

I, 


« exp[-Og(l-p(Aa-vAt)  IJ  C50) 

This  function  lias  a value  of  one  for  Aa+vAt  equal  to  zero 

2 

and  decreases  to  exp(-crg)  as  Aa+vAt  approaches  Infinity. 

It  can  now  be  seen  that  the  expI-J2(kz-^) ] phase  term  of 
Eq.  (44)  Is  cancelled  because  of  the  conjugation  of  the 
current  In  Eq.  (49).  It  also  obviously  has  no  effect  on 
the  moments  that  equal  zero. 

It  Is  necessary  to  determine  a function  to  describe 
the  normalized  covariance  function,  p(A«),  with  At  set 
equal  to  zero  for  now.  Beckman  (Ref  13:81)  suggests  a 
Gaussian  function  of  the  form 

p(Aa)  - expI-(-^-2-)2]  (51) 

c 

where  is  the  correlation  distance  defined  at  p(Ao)'»e''^ . 

For  optically  rough  surfaces  Goodman  (Ref  6:1698)  points  out 
that  Is  generally  less  than  0.1  mm.  Kurtz  (Ref  16:984) 
has  shown  that  when  of  Eq.  (29)  is  equal  to  five  or 
greater  the  correlation  distance  of  the  field  at  the  e"^ 
point  Is  then  less  than  the  correlation  distance  of  the 
rough  surface  by  a factor  of  Since  It  has  been  assumed 

that  of  Eq.  (29)  Is  greater  than  a wavelength,  then  Og 
will  be  greater  than  five  and  the  correlation  distance  of 
the  field  through  Eq.  (50)  will  be  less  than  p .02mm. 

Equation  (50)  will  now  be  called  the  field  correlation  pulse 
function  and  denoted  as  Pj^(Ao)  where  l represents  the  field 
correlation  distance.  P^(Aa)  Is  plotted  In  Fig.  7 as  a 
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function  of  ^ for  several  values  of  The  plots  clearly 
c 

show  that  the  field  correlation  distance  ts  l2S5  than  the 
surface  correlation  distance  for  the  case  of  a rough  surface 
and  a Gaussian  surface  correlation  function.  However,  this 
result  Is  not  strongly  dependent  on  the  shape  of  pCAo).  Eq. 
(50)  shows  that  any  narrow  surface  correlation  function  will 
cause  the  field  correlation  distance,  t,  to  be  very  small 
whenever  > X. 

2 2 

It  can  be  seen  by  Fig.  8 that  a -a'  Is  greatest  when 
a-a'  Is  at  Its  maximum  value,  a-a'  1s  now  limited  by  the 


O 0 

Fig.  8 Comparison  of  to  a-o' 

field  correlation  distance,  t,  and  a and  a'  maximum  Is  deter 
mined  by  the  beam  spot  size  or  the  transform  of  the  detector 
aperture,  whichever  Is  less,  through  the  terms;  hCa)  and 
h(a').  By  using  the  beam  spot  size  wCz)  as  the  limiting 
term  In  h(a) 
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Ca^-a'^)  = w^Cz)-(wCz)-i)^pwC2)j^-J^^=w  (z)j^  (52) 

'2  2 

Now  the  phase  term,  expIJ'*^^-  , can  be  approxi- 

mated by  one  If; 

z kwCz)Jl  (53) 

For  a beam  divergences,  of  less  than  three  mllllra- 

dlans  Eq,  (53)  will  hold.  Typically  beam  divergence  Is 
about  1 mllliradian,  so  the  approximation  given  by  Eq.  (53) 
appears  reasonable.  Equation  (49)  now  becomes 

00 

El  1 (t) i*{t ' )]  " //  h(a)h(ot' )a(a+vt)a(a+vt)P|,(Aa+vAt) 

-00  * 

dada'  (54) 

Equation  (54)  is  identical  to  the  far  field  case  because 
the  quadratic  terms  no  longer  affect  the  result.  At  t=t' 

Eq.  (54)  is  the  variance  of  the  complex  current.  The  cor- 
relation function  of  Eq.  (54)  is  nonstationary  because  of 
the  reflectance  terms,  so  the  correlation  distance  is  diffi- 
cult to  define.  But,  an  indication  of  the  change  in  the 
correlation  function  as  At  Is  increased  and  can  be  seen  by 
writing  Eq.  (54)  as 

00 

EI1(t)1*(t' )]  / h(a)a(o+vt)Ih(a)a(a+vt' )*p, 

,00  ^ 

(a+vAt))da  (55) 


where  * denotes  the  convolution  process. 

A graphical  picture  of  Eq,  (55)  Is  shown  In  Fig,  9,  The 
maximum  width  L of  h(a)a(a+vt)  and  h(o)a(oi+vt ' ) Is  determined 
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by  the  beam  spot  size  or  the  Fourier  transform  of  the  detec- 
tor aperture,  whichever  Is  less.  From  Eq.  (55)  and  the 
figure,  the  product  of  h(a}aCa+vtl  times  the  convolution  pro- 
cess Is  zero  once  vAt+t  Is  greater  than  L.  Because  the 
function  Is  nonstationary,  this  represents  a maximum  cor- 
relation distance  and  since  l <<  L the  correlation  distance 
Is  defined  as; 

vAt  ^ L C56) 

The  above  result  Is  rather  pleasing  since  it  would  seem 
natural  that  the  output  would  become  uncorrelated  once  all 
of  the  original  scattering  areas  passed  out  of  the  system. 
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From  Eq.  (56)  the  correlation  time  Is 

At  •«  L/v  (57) 

and  the  noise  bandwidth  Is  approximately ; 

B - 5^  " f CSS) 

The  results  of  Eqs.  (57)  and  (58)  are  similar  to  those  for 
the  case  of  filtered  "white"  noise  where  the  correlation 
time  Is  approximately  equal  to  the  response  time  of  the 
filter  and  the  noise  bandwidth  Is  equal  to  the  inverse  of 
the  filter  response  time.  However,  If  the  reflectance 
varies  much  within  the  correlation  distance  given  by  Eq. 
(56),  then  the  correlation  function  will  be  modulated  by 
the  reflectance  and  the  actual  correlation  distance  and  time 
may  become  less  than  that  given  by  the  above  equations.  The 
results  presented  so  far  have  been  based  on  the  commonly 
used  assumption  that  the  surface,  and,  thus  the  phase.  Is  a 
Gaussian  random  process.  This  assumption  has  made  the  pro- 
blem mathematically  tractible  but  Is  not  particularly 
necessary.  It  has  been  shown  that  the  characteristic  func- 
tion of  several  distributions  Is  very  small  when  Is  on 
the  order  of  a wavelength  or  greater  (Ref  14:781).  Thus, 
the  mean  of  the  complex  current  can  still  be  considered 
equal  to  zero.  Also,  the  shape  of  P^(A«)  determined  as  a 
result  of  the  Gaussian  assumption  Is  not  critical  so  long 
as  the  correlation  distance  Is  sttll  small.  This  small 
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field  correlation  distance  1s  a result  of  the  reflection  of 
a wave  from  a rough  surface,  The  rough  surface  by  defini- 
tion Is  characterized  by  a short  correlation  distance.  Thus 
the  results  presented  so  far  appear  to  apply  to  any  realis- 
tic natural  rough  surface, 

The  Signal  Detection  Models 

Since  It  has  been  shown  from  Eq.  (46)  that  the  mean  of 
the  complex  current  Is  zero,  It  follows  that  the  means  of 
the  real  and  Imaginary  parts  of  the  current  are  also  zero. 
Therefore,  the  mean  of  the  heterodyne  detector  output  cur- 
rent given  by  Eq.  (43)  Is  zero.  This  does  not  Imply  that 
the  mean  of  the  amplitude  of  the  current  Is  zero.  The 
amplitude,  A(t),  can  be  Identified  by  writing  the  complex 
current  as  the  sum  of  the  real  and  Imaginary  parts  of  the 
current  as  follows: 

1(t)  = I^+JI^.=A(t)cose(t)+JA(t)s1n0(t)  (59) 

or  by  writing  It  In  polar  from  as; 

1(t)  = A(t)e’^®^^^  (60) 

The  value  of  the  amplitude  can  be  physically  determined  from 
the  detector  output  current  by  either  coherent  signal  detec- 
tion or  envelope  detection  methods,  and  the  amplitude 
2 

squared,  A , can  be  determined  by  square  law  detection 
methods . 


L 


1; 

J 
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The  Quadrature  Model . The  statistical  models  for  de- 
tectlon  of  A and  A can  be  related  to  the  heterodyne  detec- 
tor output  current  by  first  looking  at  the  well  known 
quadrature  model  (Ref  17:238).  As  can  be  seen  from  Fig. 

10,  one  quadrature  output  Is  the  real  part  of  the  complex 
current  and  the  other  Is  the  Imaginary  part  of  the  complex 
current.  The  statistics  of  the  quadrature  outputs  can  be 


Fig.  10  The  Quadrature  Model 

determined  as  follows.  Beckman  has  shown,  by  using  the 
classical  random  walk  problem  and  the  Central  Limit  Theorem 
that  the  real  and  Imaginary  parts  of  the  scattered  field  are 
joint  Gaussian  random  variables  (Ref  5:124).  From  this  re- 
sult It  is  stral ghtforwarded  to  show  that  the  quadrature 
o,utputs  are  Gaussian  random  processes  (Ref  10:504-509). 

The  results  of  Eq.  (47)  can  be  used  to  show  that  the 
covariance  of  and  are  equal  as  follows: 


(61) 


E[i(t)1(f  )1  = ) 

+ I^(f  )I^(t)}]  = 0 
which  implies  that  for  all  t,  t' 

E[Ir(t)I^(f  )]  = E[I^(t)I^(f  )]  (62) 

and 

E[I^(t)I^(f)]  = -E[I^(f)I^(t)]  (63) 

The  terms  of  Eqs.  (62)  and  (63)  can  be  related  to  the  covar- 
iance of  the  complex  current  as  follows: 

E[i(t)i*(t')]  = E[I^(t)I^(f)  + I^.(t)I.(t')+j{I^(t)I.(t')- 

I^(f  )I.(t)}] 

= 2E[I^(t)I^(f)]+2jE[I^(t)I.(f)]  (64) 

Since  all  the  terms  in  the  expression  for  the  covariance  of 
the  complex  current  given  by  Eq.  (54)  are  real,  the  cross 
correlation  of  and  in  Eq.  (64)  is  zero  and  Eq.  (64) 
becomes 

E[i(t)i*(f)]  = 2E[I^(t)I^(f)]  (65) 

or 

E[I^(t)I^(f  ) = |E[i(t)i*(t'  )] 

2 

A of 


1 


where  k^(t,t)-l.  From  Eq.  C63l  It  can  be  seen  that  1^,  and 
are  Identically  distributed  random  processes,  In  parti- 
cular, at  t^t'  the  variances  of  the  quadratures  are  equal 
from  Eq.  (63),  and  they  are  equal  to  one  half  the  variance 
of  the  complex  current  as  shown  by  Eq.  (66).  The  results 
of  Eqs.  (64)  and  (65)  have  shown  that  and  are  uncor- 
related random  processes,  and  because  they  are  Gaussian, 
they  are  also  statistically  Independent.  Thus  at  this 
point  it  has  been  shown  that  the  quadrature  outputs  are  zero 
mean,  statistically  independent,  Identically  distributed 
Gaussian  random  processes.  Also,  the  distribution  of 
these  random  processes  has  been  related  to  the  complex 
current  by  Eqs.  (65)  and  (66). 

2 

Background  Theory  for  Detection  of  A and  A . The  above 
model  Is  now  equivalent  to  the  ^ell  known  narrowband  noise 
model  (Ref  18:399-403).  The  amplitude  squared  of  the  current 
can  be  obtained  from  the  quadrature  model  by  squaring  each 
quadrature  output  and  adding  together  the  resuits.  The 
square  root  of  this  output  Is  then  the  amplitude.  This 
model  Is  shown  in  Fig.  11,  The  first  order,  i.e.  single 
sample,  probability  density  of  the  amplitude  is  the  Rayleigh 
probability  density  which  Is  a well  known  transfori.’atlon 
from  the  Gaussian  probability  density  of  the  quadratures. 
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A(t)cos0(t)  I ^ j2 


, (.)l/2  I A(t) 


A(t)sin0(t)  I ^^^2 


Fig.  11  Determination  of  A and  A'’  from  the  Quadrature 
Outputs 


The  moment  of  the  amplitude  from  the  Rayleigh  den- 
sity Is  defined  as 

t 1) 

• rc|  t 1)  C67) 

where  r(*)  represent  the  Gamma  function.  These  Rayleigh 
moments  have  been  previously  suggested  for  a heterodyne 
system  but  the  details  of  the  model  were  not  presented 
(Ref  19:648-649). 

A typical  performance  measure  for  a system  is  the  sig- 
nal to  noise  ratio,  SNR,  The  SNR  is  defined  as  the  ratio 
of  the  signal  power  to  the  noise  power,  or,  in  statistical 
terms.  It  is  the  square  of  the  mean  divided  by  the  variance 


A SNR  of  less  than  one  implies  that  the  signal  is  dominated 
by  the  noise.  A check  of  the  SNRs  for  the  two  detection 
processes  can  easily  be  made  using  Eq.  (67)  as  follows. 


As  can  be  seen  from  the  above  two  equations  the  SNRs  are 
independent  of  all  system  parameters  and  thus  are  not  a use- 
ful performance  measure  for  system  design  purposes.  In  fact, 
the  SNR  represents  a time  independent  measure  of  performance, 
whereas  a more  useful  approach  is  to  use  the  covariance 
function  which  includes  the  time  dependence.  The  covariance 
of  A is  given  as  (Ref  18:403) 

CA(t.t')  = {2E[kQ(t.f)]-[l-k^(t.t')] 

K[k^,(t.f)]  - |}  (70) 

where  K and  E are  the  elliptical  integrals  of  the  first  and 

2 

second  kind  respectively.  The  covariance  A is  much 
simpler  and  is  (Ref  18:403) 

C^2(t,f)  = ajkj(t,f)  = [C^(t.f)]^  (71) 


w 
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The  Information  to  develop  the  individual  models  for  detec- 
2 

tion  of  A and  A is  now  available  through  the  means  from 

Eq.  (67)  and  the  covariance  functions  of  Eqs.  (70  and  (71), 

2 2 

The  Model  for  Detection  of  A . The  model  for  A will 

be  discussed  first  because  it  is  simpler  and  easier  to 

2 

Interpret.  From  Eqs.  (54)  and  (67)  the  mean  of  A can  be 
expressed  in  terms  of  the  reflectance  and  system  function 
as  follows. 

EIA^(t))=a?  = B^//h(a)h(a'  )a(a+vt)a(a'  + vt'  )Pj^(Aa)dada'  (72) 

As  previously  noted,  Pj^(Aa)  is  very  narrow  with  respect  to 
changes  in  h(a)  and  a(a)  so  as  an  approximation  it  can  be 
represented  in  terms  of  the  Dirac  delta  function  as  follows. 

Pj^(Aa)  = C6(Aa)  (73) 

60 

where  C = / Pjj^(Aa)dAa.  The  system  function,  h(a),  is  also 
-00 

symmetric  so  Eq.  (69)  becomes: 

00 

EIA^(t))  = CB^  h (-«)  h (-a)  a Ca+vt ) a (a  ' +vt  * ) <5  (a-a ' )dada ' 

OQ 

« CB^  / h^(-a)a^(atvt)d<x 
•>00 
00 

= CB^  / h^(vt-x)a^(x)dx  (74) 

•kOO 

2 

Thus  the  mean  of  A (t)  can  be  represented  by  a linear  system 
2 2 

model  of  a (x)  convolved  with  h (x),  Again  by  using  Eq,  (54), 

the  covariance  of  A from  Eq.  (71)  can  be  expressed  in  terms 

of  the  covariance  of  the  complex  current  as  follows:  ' 

j 

t 
f 
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C^2{t,t')  = [.B‘  //  h(a)hCa' )a(a+vt)a(a'+vt' )Pj^(Aa+At)dada’ ] 


[B^  //  h(vt-x)h(vt' -x' )a(x)a(x' )P, (x-x' ) 
• 00 

dxdx' 


(75) 


the  double  Integral  in  the  above  equation  is  identical  to 
that  obtained  in  computing  the  output  correlation  of  a fil- 
ter with  an  impulse  response  h(x)  driven  by  zero  mean  noise 
of  correlation  R^(x,x' ) = a(x)a(x' )P£(x-x' ) . The  square  of 
the  double  integral  in  Eq.  (75)  simply  means  that  the  noise 
can  be  represented  by  the  product  of  two  noise  processes 
that  are  identically  distributed  and  statistically  indepen- 
dent. In  contrast  to  many  other  noise  models  the  above 
model  is  signal  (i.e.  a(x) ).  dependent.  The  second  moment 
model  for  square  law  detection  is  shown  in  Fig.  12. 


A^(x) 


Bh(x) 


Bh(x) 


E[nj(x)]  = ElngCx)]  = 0 

(x.x')  = R (x,x’)  = 
n^  "2 

a(x)a(x')P-(x-x') 


Amplitude  Squared 


1 

The  speckle  noise  model  presented  above  can  be  easily 
modified  to  Include  the  quantum  noise  effects  that  were 
neglected  In  the  detector  output  current  given  by  Eq.  (21), 

It  has  been  shown  that  the  output  noise  statistics  for  a 
heterodyne  detector  can  be  represented  as  a Gaussian  random 
process  (Ref  12:189).  Thus  these  statistics  are  compatible 
with  the  model  presented  here  and  the  detector  noise  terms 
can  be  Included  simply  by  adding  the  appropriate  terms  to 
the  speckle  noise  representation  in  Eq.  (75). 

2 

Several  things  can  be  noted  from  the  model  for  A • The 
output,  A , can  be  considered  either  a spatial  signal  In 
the  variable  x or  a temporal  signal  In  the  variable  t,  where 
the  relationship  between  the  two  variables  Is  x=vt.  The 
temporal  Fourier  transform  of  A (t)  using  Eq.  (74)  Is 

FtlE[A^(t)l]  = H'(i)A'(i)  (76) 

where  H' (f )=FIh^(t) ] and  A ' ( f )=F [a^ ( t ) ] . From  Eq.  (76)  Is 

‘ can  be  seen  that  the  frequency  content  of  the  output  spectrum 

will  broaden  as  v Increases.  Thus  the  bandwidth  of  the 

electronics  that  process  the  output  must  be  scaled  by  v to 

2 

maintain  the  desired  frequency  content  a In  the  final  image. 

The  model  for  EIA  (x)]  Is  also  Identical  to  the  result 
discussed  by  Goodman  for  an  incoherent  Imaging  system  (Ref 
•7:109).  (In  two  dimensions  the  Fourier  transform  of  h (x,y) 
divided  by  the  same  transform  evaluated  at  frequencies  f^" 
f “0  Is  commonly  called  the  optical  transfer  function,  OTF.) 

y 
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j Thus,  square  law  detection  can  be  thought  of  as  an  Incoherent 

i Imaging  system  with  an  additive  noise  term  given  by  Eq.  (75). 

The  spatial  filtering  by  h (x)  and  thus  the  system 
function,  h(x),  can  be  seen  by  taking  the  spatial  Fourier 
I transform  of  E[A  (x)]  from  Eq.  (71)  as  follows: 

j Fj^tE[A^(x)]]  = B^CA'(fj^)H'(fj^) 

* B^C[A'(f^)][H(f^)*H(f^)]  (77) 

where  H(f^)=F^[h(x) ] . Since  the  Fourier  transform  has  an 
inverse  effect  on  the  width  of  a function,  the  wider  the 
system  function  in  space  the  narrower  H(fj^)  becomes.  As 

H(fj^)  becomes  narrower  so  does  the  convolution  of  H(f^)  with 

I 2 

itself  and  thus  more  of  the  high  frequency  content  of  a (x) 

is  filtered  out.  Another  way  to  see  this  filtering  effect 
is  to  look  at  the  resolution  of  two  point  sources  of  reflec- 
tivity that  are  separated  by  a distance  d,  i.e.  a^(x)=6(x)+  ^ 

6(x-d).  Then  from  Eq.  (74): 

00 

E(A^(t)]  = B^C  / h^(vt-x) [6(x)  + 5(x-d)  ]dx 
- 00 

•=  B^Ch^(vt)  + B^Ch^(vt+d)  (78) 

2 

As  an  example,  it  can  be  assumed  that  h (x)  is  a Gaussian 
amplitude  function  (consistent  with  a Gaussian  laser  beam), 
and  then  the  resolution  limit  can  be  picked  as  the  point 
where  the  peak  of  one  function  is  at  the  e"^  point  of  the  1 

I 

other.  Then  from  Eq.  (78)  it  can  easily  be  seen  that  d must 
be  the  distance  at  which  h^(x)  is  at  it's  e"^  point  if  the 
two  points  are  to  be  resolved.  So  again  the  wider  system  j 
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function  in  space  the  worse  the  resolution  ability  of  the 
system.  This  is  a very  Intuitive  result  since  the  detector 
measures  an  average  of  the  reflectance  points  each  weighted 
by  the  value  of  the  system  function  at  that  point.  The 
wider  h(x)  1s  then  the  more  even  this  weighting  becomes  over 
the  width  of  h(x).  Thus,  If  two  points  of  equal  reflec- 
tivity are  to  be  resolved  by  the  system,  they  must  be 
separated  by  a distance  large  enough  such  that  the  dif- 
ference In  the  weighting  of  each  point  Is  large  enough  to 
cause  a distinguishable  effect  In  the  detector  output. 

The  effect  of  h(x)  on  the  mean  does  not  tell  the  com- 
plete  story  since  h(x)  also  effects  the  covariance  of  A , 
le.e  the  fluctuations  of  A Cx).  The  covariance  of  A from 
Eq.  (75)  Is  nonstationary,  but  an  indication  of  the  correla- 
tion distance  can  be  determined  from  the  correlation  distance 
of  the  complex  current  given  by  Eqs.  (55)  and  (56).  The 
result  of  squaring  the  complex  current  correlation  function 
Is  to  reduce  the  maximum  correlation  distance  of  Eq.  (56) 
slightly.  Still,  It  can  be  seen  that  the  wider  the  system 
function  the  longer  the  maximum  correlation  distance. 

The  square  root  of  the  covariance  at  t-t  from  Eq.  (75), 
I.e.  the  square  root  of  the  variance,  represents  the  rms 
variation  from  the  mean  due  to  the  "noise".  This  variation 
can  the  thought  of  as  the  contrast  variations  In  the  output 
Image  that  are  due  to  the  rough  surface  (speckle)  noise. 

Also,  the  correlation  distance  represents  the  average  period 
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over  which  this  noise  process  is  related  or  does  not  change 
"very  much".  Thus,  the  the  correlation  distance  can  be 
thought  of  as  the  average  speckle  cell  size  in  the  resulting 
image,  and  it  is  directly  related  to  the  width  of  the  system 
function,  h(x),  as  was  shown  by  Eq.  (56). 

The  effect  of  the  system  function  on  the  noise  can 
also  be  seen  by  looking  at  the  noise  power  spectral  density. 
The  Fourier  transform  of  the  output  covariance  with  respect 
to  vAt  is  called  the  power  spectral  density,  S(fj^)  (Ref  20: 
347).  It  represents  the  noise  power  per  spatial  frequency 
that  passes  through  the  system.  The  total  output  noise 
power  is  found  by  integrating  S(f^)  over  all  frequencies. 

The  power  in  a particular  frequency  band  is  found  by  simply 
integrating  S(f^)  over  that  band  of  frequencies. 

The  covariance  of  A given  by  Eq.  (75)  is  nonstationary 
so  strictly  speaking  it  is  not  subject  to  Fourier  analysis. 
But  if  the  reflectance  is  considered  nearly  constant  (or 
slowly  varying)  then  the  power  spectral  density  of  A is 

S 2(fj^)  = B^a^Pj^p(f^)|H(f J|2*B^a^P^p(f^)|H(fJ|2  (79) 

A 

where  P^p(fj^)  = FtPj^(Ax)  ] . Because  Pj^(Ax)  is  very  narrow  com- 
pared to  h(x),  the  spatial  frequency  spectrum  of  P£p(^j() 
will  be  wideband  or  "white"  compared  to  H(f^).  Or,  by  again 
approximating  Pjj^(Ax)  as  C6(Ax)  as  in  Eq.  (73),  Eq.  (79)  can 
be  written  as: 
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(80) 


2 

Thus,  as  h(x)’  becomes  wider,  lH(fj^)|  becomes  narrower,  and 

the  total  noise  power  out  will  decrease  because  of  the 

smaller  non-zero  ranges  over  which  S o(f„)  will  exist. 

* 2 

This  completes  the  analysis  of  the  A model.  It  has 
been  shown  by  several  approaches  that  the  resolution  of  the 
reflectance  information  is  reduced  as  the  system  function  is 
widened  to  reduce  the  noise  power  or  noise  flucuations. 

The  model  for  envelope  detection  will  now  be  developed  to 
see  how  the  system  function  affects  that  method  of  detection. 

The  Model  for  Detection  of  A.  The  model  for  the 
measurement  of  the  amplitude  is  much  more  difficult  because 
of  the  expression  for  the  covariance  and  the  square  root 
involved.  The  mean  of  the  amplitude  from  Eqs.  (54)  and 
967)  is: 


E[A(t)]  = B-|  [//h(a)h(a' )a(a+vt)a(a'+vt)Pj^(Aa)dada' 

» C^'^^[/h^{a)a^(a+vt)dx] 

= B_n  jgjj 

This  is  the  same  expression  (within  a constant)  as  the  square 
root  of  the  mean  of  A from  Eq.  (74).  Therefore  the  comments 
following  Eqs.  (76)  and  (77)  regarding  the  electronics  band- 
width and  the  system  function  apply  for  the  detection  of 
the  amplitude,  too. 
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The  covariance  of  A from  Eq.  (70)  was  expressed  in 
terms  of  the  elliptical  integrals  of  K and  £ as  follows: 

o? 

CA(t.f)  = -2’  ^^ElkpCt.f  )]-[l-k‘(t.t' )]K[k^j(t,f  )]-f>  (70) 

The  elliptical  integrals  can  be  represented  by  the  following 
series  (Ref  21:310) 

IC(X)  = f (82) 


1.3.5. ..(2j-i), 2 x' 
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The  terms  in  the  outer  brackets  of  Eq.  (70)  can  be  combined 
using  x=k^(t,t')  as  follows: 

, 1 . 3. 5. ... (2 j- 1 ) , 2 „ 2 J 1 n 


2*4*6* • • • (2  J, 


]"  x"'}  - J 


j-1  ^ 4*6 (2j)  ^ ^^2j-r  ^2^  4j-2J^ 

Eq.  (84)  can  now  be  evaluated  explicitly  for  the  first  terms 
in  the  series  as  follows 


8 ^ ^8'16'  ^ * fi  'fia'  ^ ^ ^ 
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SO  the  covariance  of  the  amplitude  becomes: 
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The  first  term,  except  for  the  constants.  Is  the  same  as 
the  covariance  of  A from  Eq.  (71),  where  the  square  of 
C^(t,t')  implied  that  the  noise  could  be  modeled  as  the 
product  of  two  Independent  noise  processes.  In  Eq.  (86) 
the  sum  of  the  higher  order  powers  of  C^(t,t')  Implies  that 
the  noise  can  be  molded  as  a sum  of  Independent  noise  pro- 
cesses. Additionally  each  term  in  the  sum  is  Itself  a 
product  of  m Independent  noise  processes  where  m Is  equal 
to  the  power  of  the  C^(t,t')  term  Involved. 

For  any  t'  not  equal  to  t In  Eq.  (86)  the  term  C|(t,t') 

p 

Is  less  than  so  the  terms  of  Eq.  (86)  are  rapidly  de- 


creasing with  the  higher  powers  of 


C^(t.t' ) 


As  the  CO- 


I 

herence  length  of  C^(t,t')  of  Eq.  (56)  is  approached 

C.Ct.f)  ^ ^ 

~ — 2 becomes  very  small  and  Eq.  (86)  can  be  quite 

accurately  approximated  by  the  first  term.  Except  for  the 
constants  this  approximation  yields  the  same  form  as  the 
covariance  of  A from  Eq.  (75),  Thus  the  coherence  dis- 
tance can  be  considered  to  be  approximately  the  same  for 
both  methods  of  detection. 
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If  the  power  spectral  density  of  Eq.  (86)  Is  deter- 
mined  in  a manner  similar  to  the  case  for  A from  Eq,  (80), 
the  result  would  be  a series  of  double,  quadruple,  and 
higher  order  convolutions.  This  means  that  the  noise 
power  is  unlimited  in  frequency,  although  it  is  reduced  in 
amplitude  at  each  higher  order  convolution.  In  fact,  if 
the  value  of  the  first  term  of  Eq.  (86)  at  t=t'  is  com- 
pared to  the  variance  of  A using  Eq.  (67),  the  result  is: 


(Jn, 
' 16  > 


2(2-f) 


- ( ^ ) = 915 

8 M-n^ 


(87) 


Eq.  (87)  shows  that  91.5%  of  the  noise  power  is  concentrated 
in  the  first  term  of  the  covariance  of  A when  the  covariance 
function  is  at  it's  maximum  value,  i.e.  at  t=t'.  It  was 
also  determined  in  the  previous  paragraph  that  the  first 
term  was  an  excellent  approximation  for  the  covariance  of 
A when  the  covariance  is  small,  i.e.  at  the  coherence  time. 
Thus,  as  a first  approximation  the  noise  covariance  func- 
tion can  be  represented  as: 

,2 


C.(t,t')  = 


n[C^(t,t' )]' 


16  a 


(88) 


i 


for  all  t,t'.  This  approximation  is  in  addition  to  the  more 
accurate  first  term  approximation  that  was  given  for  the  co- 
herence length.  Again,  except  for  the  constants,  the  npise 
representation  is  now  the  same  as  that  for  A given  by  Eq. 
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(71),  With  this  approximation  the  model  for  detection  of 
A Is  given  in  Fig.  13. 


Fig.  13  The  Approximated  Model  for  Detection  of 


the  Current  Amplitude 

The  system  parameters  now  have  the  same  effect  on  the 

resulting  Image  in  both  the  A and  A models  although  the 

final  image, is  not  the  same  because  of  the  square  root  of 

the  mean  In  the  model  for  detection  of  A.  In  fact,  the 
2 2 

mean  of  A Is  proportional  to  a , I.e.,  the  square  of  the 
field  reflectivity,  and  thus  Is  directly  related  to  the  in- 
tensity of  the  reflected  field.  Since  the  human  eye  also 
measures  Intensity,  the  A Image  Is  proportional  to  the 
reflectivity  term  that  people  are  familiar  with.  On  the 
'Other  hand,  the  mean  of  A Is  proportional  to  the  field  re- 
flectivity, a,  and  the  result  is  an  image  that  may  not  be 
natural  to  the  human  observer.  However,  this  image  is  still 
a valid  representation  of  the  surface  characteristics. 
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IV.  SUMMARY 


Conclusions  ~ 

As  was  stated  in  the  background  section,  there  is  a 
great  deal  of  information  available  in  the  literature  con- 
cerning speckle  as  it  relates  to  imaging  systems  that  make 
intensity  measurements.  Very  little  information  was  avail- 
able on  modeling  the  output  of  an  optical  heterodyne  line- 
scan  imaging  system.  The  purpose  of  this  thesis  was  to 
develop  such  a model  that  included  the  effects  of  the 
speckle  noise,  which  is  caused  by  reflection  of  the  field 
from  a rough  surface  object.  This  rough  surface,  however, 
is  necessary  in  the  scanning  system  to  provide  backscat- 
tered  radiation  to  the  detector,  otherwise  the  detector 
would  only  receive  a return  signal  the  few  times  the  laser 
beam  was  normally  incident  to  the  object  surface.  The  met- 

■3 

¥ 

hods  of  detecting  either  the  amplitude  or  amplitude  squared  ] 

of  the  current  from  the  heterodyne  detector  have  been 
modeled  in  this  thesis.  The  amplitude  of  the  current  is 
easily  measured  in  an  actual  system  by  envelope  detection 
or  it  can  be  detected  by  more  sophisticated  coherent  detec- 

i: 

tipn  methods.  The  amplitude  squared  of  the  current  can  be 

determined  by  use  of  common  square  law  detectors.  | 

I 

Because  the  detector  output  is  the  average  of  all  of  ji 

the  reflectance  points  within  the  laser  beam,  it  was  argued 
that  the  significant  changes  in  the  detector  output  occurred 
as  a function  of  the  scanning  direction  of  the  system.  The 
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reflectance  at  each  point  in  the  scanning  direction  could 
be  considered  equal  to  the  average  of  the  reflectance 
across  the  width  of  the  beam  perpendicular  to  the  scanning 
direction  at  that  point.  Thus,  the  system  model  was  devel- 
oped in  the  scanning  direction  only,  The  one  dimensional 
Huygens-Fresnel  integral  was  used  to  propagate  the  laser 
output  to  the  object's  surface  and  the  reflected  field  back 
to  the  detector.  The  laser  and  detector  were  modeled  as 
being  stationary  in  a fixed  coordinate  system  while  the  sur- 
face  characteristics  were  moved  beneath  the  laser  beam.  It 
was  pointed  out  that  the  effects  of  beam  spreading  and 
changing  surface  velocity,  that  were  ignored  in  this  model 
but  which  occur  in  the  actual  system,  could  be  determined 
from  this  model  simply  by  varying  the  appropriate  parameters 
In  the  development  of  the  system  model,  the  object's  surface 
was  considered  to  be  a zero  mean  Gaussian  random  process, 
but  it  was  later  argued  that  the  results  were  typical  of 
any  naturally  rough  surface.  It  was  shown  that  the  rms 
roughness  of  the  surface  could  be  compared  directly  to  the 
optical  wavelength  and  that  many  surfaces  are  rough  com- 
pared to  optical  wavelengths.  A system  function  was  defined 
which  represented  the  combined  effect  of  the  detector  field 
of  view  and  the  laser  beam  amplitude  function,  Then  the  com 
plex  current  was  used  in  determining  the  statistics  of  the 
optical  detector's  output  and  it  was  shown  that  the  mean 
signal  was  zero  when  the  objects  surface  was  rough  compared 


to  an  optical  wavelength.  It  was  also  shown  that  the  cor- 
relation function  of  the  reflected  field  was  very  narrow, 
on  the  order  of  tenths  and  hundredths  of  a millimeter.  This 
narrow  correlation  function  meant  that  the  reflected  field 
was  now  spatially  coherent  over  just  a very  short  distance. 
This  fact  allowed  the  quadratic  phase  terms  to  be  neglected 
which  resulted  in  a model  that  was  the  same  for  both  far 
field  and  near  field  cases. 

Because  the  mean  of  the  detector  output  current  was 
zero,  a mehtod  of  detection  the  amplitude  or  amplitude 
squared  of  the  current  was  included.  The  amplitude  is  re- 
lated to  the  sum  of  the  squares  of  the  real  and  imaginary 
parts  of  the  current  so  the  quadrature  model  was  presented. 
The  quadrature  outputs  were  shown  to  be  identically  distri- 
buted, zero  mean  joint  Gaussian  random  processes.  Based  on 
this  result,  it  was  concluded  that  the  first  order  density 
of  the  amplitude  was  Rayleigh  and  that  the  model  was  now 
similar  to  the  well  known  narrowband  noise  model.  The  sig- 
nal to  noise  ratio  for  each  detection  method  was  calculated 
from  the  moments  of  the  Rayleigh  density  but  in  both  cases 
it  was  Independent  of  the  system  parameters. 

However,  a second  moment  model  for  each  detection  pro- 
cess, which  included  the  system  parameters,  was  developed. 

The  respective  mean  (signal)  and  covariance  (noise)  functions 
were  all  expressed  in  terms  of  the  previously  determined  co- 
variance  function  of  the  complex  current.  The  amplitude 
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squared  signal  model  was  shown  to  be  identical  to  the  linear 
system  model  for  an  Incoherent  imaging  system  where  the 
Fourier  Transform  of  the  system  function  squared  represents 
the  well  known  Optical  Transfer  Function.  It  was  shown 
that  the  resolution  ability  of  the  system  degraded  as  the 
width  of  system  function  was  increased  in  space.  This  was 
Intuitively  pleasing  since  an  increase  in  the  width  of  the 
system  function  causes  more  reflectance  points  to  be  in- 
cluded in  the  averaging  process  and  thus  the  effect  of  each 
individual  reflectance  point  becomes  smaller.  The  amplitude 
signal  model  was  developed  and  shown  to  be  identical  to  the 
amplitude  squared  signal  model  except  for  some  constants 
and  the  final  square  root  of  the  output.  The  filtering 
effects  of  the  system  function  were  the  same  for  both  models. 

In  each  model,  the  noise  covariance  was  shown  to  be 
represented  by  the  square  of  the  complex  current  correlation 
function  and  unlike  some  common  noise  processes  it  was  a 
function  of  the  reflectivity  signal.  The  models  provided 
considerable  information  about  the  speckle  noise.  The  co- 
herence length  of  the  noise  process  represents  the  period 
that  noise  is  related  and  thus  corresponds  to  the  average 
speckle  cell  size  in  the  image.  This  correlation  distance 
was  determined  to  be  equal  to  the  width  of  the  system  func- 
tion except  that  when  the  reflectivity  varied  significantly 
over  the  same  distance,  it  could  become  smaller.  The  con- 
trast in  the  image  due  to  the  speckle  noise  was  given  by 
the  square  root  of  the  variance  and  could  be  determined  from 
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the  noise  representation  in  the  models.  The  models  show 
how  the  system  parameters  affect  the  signal  and  the  speckle 
noise  and  they  provide  the  basis  for  developing  signal  pro- 
cessing methods  that  will  optimize  the  desired  image.  Also, 
it  was  shown  that  the  common  heterodyne  detector  noise 
models  could  be  easily  added  to  the  speckle  noise  models 
presented  here. 

Recommendati ons 

Although  the  models  developed  in  this  thesis  describe 
the  speckle  effects,  they  do  not  provide  a definite  solution 
to  the  problem  of  producing  the  "best"  image  in  the  pre- 
sence of  speckle  noise.  It  is  recommended  that  optimization 
techniques  be  applied  to  the  models  to  determine  what  signal 
•processing  could  be  done  to  provide  satisfactory  system  per- 
formance. It  would  be  helpful  to  know  what  the  expected 
spatial  frequency  content  of  the  reflectance  is  so  that  the 
reflectance  signal  power  spectrum  could  be  compared  to  the 
noise  power  spectral  density.  Then  a minimum  acceptable 
resolution  criterion  could  be  developed,  and  an  appropriate 
filter  determined  that  would  yield  the  desired  resolution 
while  filtering  those  frequencies  that  reduce  the  image 
quality.  Also,  this  would  allow  comparisons  to  be  made  be- 
tween the  two  signal  detection  models  to  determine  which  de- 
tection method  would  yield  the  best  performance.  Of  course 
comparison  of  the  models  developed  here  with  additional 


experimental  data  would  serve  to  validate  them.  Finally, 

60 


L 


the  development  of  a similar  model  for  optical  direct  de- 
tection line-scan  Imaging  systems  would  provide  the  basis 
needed  for  determining  which  type  of  system  should  be 
further  developed. 
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